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Introduction .— Predicting properties of hadrons still 
represents a serious challenge for Quantum Chromody¬ 
namics (QCD). Heavy quarks closely resemble static test 
charges and therefore are useful to probe confining prop¬ 
erties of QCD. So far, great progress has been made in 
the study of quarkonia, i.e. mesonic states that contain 
two heavy constituent quarks. In contrast, systems of 
three or more heavy quarks, which are a good starting 
point for understanding the phenomenology of baryons 
and multi-quark bound states, are much less studied. In 
this case a key issue is whether multi-quark interactions 
can be understood in terms of two-body interactions or 
whether there are genuine three- and many-body effects 
to be considered as part of the overall picture of strong 
interactions mm- 

The best known phenomenological models of the N- 
quark potential are those of A = 3, the so called A and 
Y-laws [3]. The A-law is based on pairwise interactions 
between quarks, while the Y-law is an examples of three- 
body interactions. In the infrared region the former pre¬ 
dicts that the potential grows linearly with the perimeter 
of the triangle formed by quarks [1] , while the latter pre¬ 
dicts a linear growth with the minimal length of a string 
network which has a junction at the Fermat point of the 
triangle [5]. 

Until recently, lattice gauge theory was the premier 
method for obtaining quantitative and qualitative infor¬ 
mation about strongly interacting gauge theories. For 
the three-quark potential the accuracy of numerical sim¬ 
ulations has been improved during the past decade Em 
that provided evidence for the Y-law at long distances. 
On the other hand, it is expected that at short distances 
the A-law is a good approximation to the potential mm- 
However, what is still missing is a model which would in¬ 
corporate the A-law at short distances and the Y-law at 
long ones. 

In this Letter we present the first example of such 
a model. It continues a series of studies [lUHH] de¬ 
voted to the static potentials in four-dimensional (pure) 
gauge theory by means of a five (ten)-dimensional effec¬ 
tive string theory. Our reasons for continuing to pursue 
this model are: 

(1) Because there is no string theory which is dual to 


QCD. It would seem very good to gain what experience 
we can by solving any problems that can be solved within 
the effective string model already at our disposal. 

(2) Because the results provided by this model are 
consistent with the lattice calculations and QCD phe¬ 
nomenology [SHin]. 

(3) Because analytic formulas are obtained by solving 
this model. 

(4) Because it allows us to make predictions [16] which 
may then be tested by means of other methods, e.g., nu¬ 
merical simulations. 

Before proceeding to the detailed analysis, let us set 
the basic framework. As for the quark-antiquark po¬ 
tential, the static A-quark potential can be determined 
from the expectation value of a Wilson loop. The loop 
in question, baryonic loop, is defined in a gauge-invariant 
manner as IYnq = with the 

path-ordered exponents along the lines shown in 

Figure 1. In the limit T —> oo the expectation value of 
the loop is simply (IYnq(C')) ~ e“'®'^, with E the ground 
state energy of A quarks (A-quark potential). 

In discussing baryonic Wilson loops, we adapt the for¬ 
malism [niiiH] proposed within the AdS/CFT correspon¬ 
dence m to our purposes. First, we take the following 
ansatz for the background geometry |5D| 

= e^’' ^{dt^ -P dx^ + dr'^) -P doj°'doj^ , (1) 

where dx'^ = dx"^ -k dy"^ + dz^. This is a deformed prod¬ 
uct of AdSs and an internal space (five-sphere) X whose 
coordinates are a;“. The deformation is due to the r- 
dependent warp factor, with s the deformation parame¬ 
ter. Such a deformation is a kind of the soft wall model of 
[21j . where the violation of conformal symmetry is man¬ 
ifest in the background metric. In Q, there are two free 
parameters to be fitted to the results of numerical sim¬ 
ulations or quarkonia spectra. Both fits look very good 

[13 Ej. 

Next, we consider the baryon vertex which is a A- 
string junction. Since we are interested in a static quark 
potential, we choose static gauge and then make an 
ansatz for the action, describing a static configuration, 
of the form 
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FIG. 1: Left: A baryonic Wilson loop in SU{S) gauge theory. Right: In SU{4), a configuration used to calculate the 
expectation value of a baryonic loop. The quarks are set on the x-y plane. R is a baryon vertex located at r = ro and 5 is its 
projection onto the x-y plane. 


Q-25r^ 

5vert = m -T, (2) 

r 

where m and s are parameters, r is independent of t, and 
T = dt. In what follows, we will assume that quarks 
are placed at points on the boundary of 5-dimensional 
deformed AdS (at r = 0) but at the same point in the 
internal space. This assumption makes the problem effec¬ 
tively five-dimensional. Therefore the detailed structure 
of X is not important, except for the warp factor de¬ 
pending on the radial direction. The motivation for such 
a factor in ([^ is drawn from the AdS/CFT construction, 
where the baryon vertex is a 5-brane m- Taking a term 
f dtd^uj\fg^ from the world-volume action of the brane 
results in Te“^®’'^/r if r is independent of t. This is, 
of course, a heuristic argument but, as we will see, the 
ansatz ([^ is quite successful: it allows us to describe the 
results for iV = 3 using just one parameter. 

The expectation value of the Wilson loop is schemati¬ 
cally given by the path integral over world-sheet fields 

(f^NQ(C)) = J , (3) 

where Sw is a total action of the Nambu-Goto strings 
and vertex. The strings are stretched between the quarks 


on the boundary and the baryon vertex in the interior, 
as sketched in Figure In principle, the integral can 
be evaluated approximately in terms of minimal surfaces 
that obey the boundary conditions. The result is writ¬ 
ten as (Wnq(C')) = exp[—where Sn means a 

renormalized minimal area whose weight is Wn- 

Calculating the N-quark potential. — We consider a sit¬ 
uation in which N quarks are placed at the vertices of a 
regular Wsided convex polygon of side length L. This 
configuration has the symmetry group Hence 5 is a 
center of the polygon and all the strings have an identi¬ 
cal profile. To compute the potential, we proceed along 
very similar lines to those of [12]. First, we take the static 
gauge that allows us to solve the equations of motion and 
determine the string profile. Next we extremize the ac¬ 
tion with respect to the location of the baryon vertex rp 
that results in the no-force condition at r = rg. There is, 
however, one important distinction between the present 
calculation and those in the literature devoted to large N 
gauge theories. We make an assumption that the param¬ 
eter m is negative. As a result, gravity pulls the vertex 
toward the boundary. This bends the strings and blunts 
the tip of the configuration |16j . as shown in Figure 1. 

Having found the solution, we can compute the total 
energy of the configuration. At the end of the day we 
arrive at [16] 


L{v) = 2sin(^) 


TT \ /A 


UO 


dv (l - " + r du (l - ’ 


*(> 


and 


E{y)=NQ^j- 




( 4 ) 


+ C, (5) 


where v = srg, g = 2 ^^, k = and C is a nor- malization constant. A is a function of and/ 1 : such that 
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A = —ProductLog[—z/e“'^(l —+ where 

ProductLog(z) is the principal solution for w in z = we^ 
[ 22 ] , Also note that u G [0, i^*], with i/* a solution to 

Thus, the A^-quark potential as a function of the in¬ 
terquark separation is parametrically specified via the 
parametric functions E{i') and L{v). Importantly, the 
parameters 0 and s coincide with those of the quark- 
antiquark potential mm and, as a consequence, n is 
the only free parameter in the model. This is our main 
result. 

It is worth analyzing E{L) in the two limiting cases, 
short and long distances. In the former case we find 

E{L) = — (? N —^ + C + ^ctqL + o{L) , ( 6 ) 


OiNQ — - LqEqQ, (To — -r— (Ei + —^Eo)q5, (7) 

where T/r^, Vn = r,^ denotes 

the distance between the vertices i and j, and E^ are 
given in the Appendix. 

In the latter case, we get a generalization (star-law) of 
the Y-law with a single Steiner point S 

N 

„ • / TT \ ’ (^) 

2 sm(^) 

with the same string tension a as in [Mia 


CT = egs, 



dv 




C. (9) 


Three features of the model are worth highlighting here. 
First, at short distances it yields the subleading linear 
term [^. Second, at long distances the model reduces 
to the star-law with the physical string tension ct, as ex¬ 
pected [12]. Finally, the constant terms at short and long 
distances are different. Notice that c — C is scheme in¬ 
dependent and is free from divergences. This makes the 
model so different from the phenomenological laws. 

Numerical Results and Phenomenological Prospects .— 
It is of great interest to compare our model of the N- 
quark potential with the results of numerical simulations. 

Clearly, A^ = 3 is of primary importance [24]. In this 
case, we set g = 0.176 and s = 0.44GeV^, i.e., to the 
same values as those of m used for modeling the quark- 
antiquark potentials of dSj. With these parameters fixed, 
the model has only a single free parameter k. We then fit 
it with a^Q = 0.125 taken from [ 5 ]. So we get k = —0.083. 
The result is plotted in Figure 2, on the left. We see that 
the model reproduces the lattice data remarkably well. 
Also note that a fit to Q! 3 q = 0.129 of [3] doesn’t change 
the picture as the discrepancy between the data of [ 5 ] 
and [3] is negligible at small distances. 

For completeness, it is worth making a couple of esti¬ 
mates. First, from Q we get Q! 3 q/q;qq « 0.495, where 
ttQQ = (27r)^r“^(|;)g [ini |2S]. Thus, the relation be¬ 
tween the ’’Coulomb” coefficients found in perturbative 
QCD holds with good accuracy in our model. This looks 
puzzling as we consider small distances but not very small 
ones of perturbative QCD. Second, from Q and (|^, we 
find (To/cr « 1.007 that favours the A-law at short dis¬ 
tances, as also noted in [ 8 ]. 

For practical purposes, the parametric form of the po¬ 


tential looks somewhat awkward. It is instructive to com¬ 
pare the lattice data to the asymptotic behavior of E{L) 
to see what happens. In Figure 2 we have plotted the 
results. As can be seen, in the range of interest a single 
Coulomb-type term doesn’t yield a satisfactory descrip¬ 
tion. But if one adds an additional linear term, then the 
situation will improve. Such a two-component model, al¬ 
most the A-law, does describe the data quite well in the 
range 0.1 fm < L < 0.6 fm. However, at longer distances 
0.6 fm < A < 1.2 fm, it becomes less accurate than a sin¬ 
gle linear term from the Y-law, as already noted in [7]. 
Thus the model we propose smoothly interpolates be¬ 
tween the A and Y-laws with a transition at A « 0.6 fm. 

Unfortunately, in lattice gauge theory very little is 
known about the A-quark potential if A > 3. Even in 
SU(4) [5], the availability of data is much more limited 
than it should be to make a consistent comparison with 
our predictions. 

In the physically interesting case A = 3, the model 
incorporates pairwise interactions, the A-law, at short 
distances and a genuine three-body interaction at longer 
ones. Note that one can think of the pairwise interac¬ 
tion as that of a quark and a diquark thanks to a string 
stretched in between. Such an interaction occurs at inter¬ 
mediate scales in between very small scales, where quarks 
are asymptotically free, and large scales, where quarks 
are strongly coupled. 

Summary and discussion .— In this Letter, we have 
modeled the A-quark potential using the now standard 
ideas motivated by gauge-string duality. Our work based 
on the background geometry which is singled out 
by the earlier works [muD], provides the first convinc- 
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FIG. 2: The lattice data, obtained on an equilateral triangle at /3 = 6.0, are taken from (6] [8] (squares), [7] (disks), and [9] 
(triangles). We use the normalization of [S]. We don’t display any error bars because they are comparable to the size of the 
symbols. Left: The 3-quark potential at g = 0.176, s — 0.44GeV^, k = —0.083, and C = 1.87GeV. Right: Some asymptotic 
curves for the model at the same parameter values: —-I- C (dot dashed), —3^^D^ + C + |(JoI/ (long dashes), and y/SaL + c 
(short dashes). 


ing example of interpolation between the A and y-laws. 
Mathematically, the potential is described by a compli¬ 
cated function whose asymptotic behavior is given by the 
A and M-laws. 

The model we are developing is an effective string the¬ 
ory based on the Nambu-Goto formalism in a curved 
space. Therefore it has some limitations including: the 
issue of a Liischer-like correction on a curved background 
and the issue of attraction between the baryon vertex and 
boundary. What could be the reason for the latter? m 
is a result of a resummation of infinitely many terms {a' 
corrections) in the five-brane action. Is it negative be¬ 
cause the brane tension is negative [5S], and if so, does it 
lead to instability? These questions have no obvious or 
immediate answers. Hopefully, it will be resolved in the 
future by using the Green-Schwarz formalism, already de¬ 
veloped for strings on AdSs x S® [27]. Obviously, finding 
the way to the string description of QCD is a challeng¬ 
ing and difficult problem. In the meantime, lattice gauge 
theory and effective string models will remain the main 
tools of investigation. 
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Appendix: Some useful formulas .— In studying the 
short distance behavior of E{L), the following two facts 
are useful. First, L goes to zero as zz —>• 0. Second, for 
small 12 and A we have zz(A) = p^\+ {l — 2 k^ — + 

o(A^), with p = I — 


Expanding and to subleading order in A, we 
find 


+ E — +EiXXj+C, (A.l) 


L = 
together with 


Lo = -sin(^)H(K2;i, I) , 


L, = ^sin(^)(H(.^;-i,|)-H(.^;-i,|)-,), 
E, = n(-2kp^ + i, i) + ■ 


Here p = 2pi |k|“^(1 — — P®), B{z; a, b) = B{a, b) + 

B{z; a, b), and B{z-, a, b) is the incomplete beta function. 
Then a simple algebra leads to the explicit formulas for 
Onq and CTo- 
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